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Abstract 

The emission of particles from black holes created in the early 
Universe has detectable astrophysical consequences. The most 
stringent bound on their abundance has been obtained from the 
absence of a detectable diffuse flux of 100 MeV photons. Fur¬ 
ther scrutiny of these bounds is of interest as they, for instance, 
rule out primordial black holes as a dark matter candidate. We 
here point out that these bounds can, in principle, be improved 
by studying the diffuse cosmic neutrino flux. Measurements of 
near-vertical atmospheric neutrino fluxes in a region of low geo¬ 
magnetic latitude can provide a competitive bound. The most fa¬ 
vorable energy to detect a possible diffuse flux of primordial black 
hole origin is found to be a few MeV. We also show that measure¬ 
ments of the diffuse u T flux is the most promising to improve the 
existing bounds deduced from gamma-ray measurements. Neu¬ 
trinos from individual black hole explosions can be detected in 
the GeV-TeV energy region. We find that the kilometer-scale 
detectors, recently proposed, are able to establish competitive 
bounds. 


1 Introduction 


Soon after Hawking described the quantum evaporation of black holes [p]], 
it was suggested that emission from black holes, created in the early Uni¬ 
verse, has detectable astrophysical consequences |2], ffj. Such black holes 
are referred to as primordial. They may be formed by the collapse of den¬ 
sity fluctuations at an early epoch [|], at a cosmic phase transition or 
by more exotic mechanisms such as the collapse of closed cosmic strings [j7j]. 
Attempts were made to search the Hawking radiation from primordial black 
holes (PBHs). The instantaneous emission spectra [|lj of particles by an indi¬ 
vidual black hole resemble a black body spectrum with temperature inversely 
proportional to the black hole mass. Only holes of masses below the critical 
mass M* ~ 4 — 6 x 10 14 g evaporate within the age of the Universe through 
this mechanism. 

In this paper we assume that a particle species is emitted when the black 
hole temperature exceeds the relevant rest mass and if the particle appears 
elementary at this temperature. A solar mass black hole emits only massless 
particles, photons and possibly neutrinos. As a black hole evaporates its 
temperature rises crossing the emission threshold of heavier particles so that 
new degrees of freedom are added to the emission. A black hole of mass M* 
also emits electrons and positrons with energies in the 50 — 100 MeV energy 
range. This average energy is sufficiently close to the muon and pion rest 
mass so that they are also kinematically accessible. In the Standard Model 
this means that the quark and gluon degrees of freedom will contribute to 
the radiation. When the temperature approaches A q C d quarks and gluons 
are emitted freely which subsequently fragment into hadrons by mechanisms 
familiar from accelerator studies. Fragmentation products decay into stable 
particles which add to the direct Hawking emission, thus contributing an 
enhancement in the spectrum near ~ 140 MeV, the pion rest mass [|5|. 

PBHs can be searched for as nearby point sources or, alternatively, one 
may search for the accumulated diffuse flux of all PBHs which evaporated 
within the lifetime of our Universe. Page and Hawking obtained the first 
such bound from the requirement that this diffuse flux cannot exceed the 
observed gamma ray spectrum |J. The most stringent bound is obtained 
for photons with energy of order 100 MeV. They only considered the direct 
photon flux. This result was later updated to include all Standard Model 
degrees of freedom as well as the fragmentation of quarks and gluons |J. 
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Bounds obtained from the observed diffuse electron, positron and antiproton 
spectrum [fjj are remarkably close to those from gamma rays. 

Searches for individual, nearby holes provide bounds on explosion rates 
which are complementary to the diffuse bounds. These bounds depend on 
the possible clustering of PBffs around galaxies and on their final stage 
evolution which is dictated by yet unsettled high energy physics. Though 
the solutions to these problems are somewhat ambiguous, an observation of 
such an explosion would reveal all the degrees of freedom of all particles in 
the Universe irrespective of their mass. 

The extension of the work on photons to include neutrinos is most nat¬ 
ural. The construction of a new generation of solar and supernova neutrino 
detectors as well as the commissioning of a new type of high energy neutrino 
telescope provide us with novel ways to search for PBHs. The low cross sec¬ 
tions for neutrino interactions require large detector volumes, but they also 
allow detection on the Earth’s surface avoiding shielding by either the inter¬ 
stellar medium, the atmosphere or the Earth itself. The emission of direct 
neutrinos has been previously addressed by Carr [Q]. In this work we include 
the production of neutrinos via the decay of hadrons. By comparing the 
diffuse neutrino fluxes from PBffs to the observed neutrino backgrounds we 
determine the most favorable energy regions for obtaining new independent 
bounds on PBff abundances. We will show that such regions of opportu¬ 
nity exist in the range of sensitivity of both low- and high energy neutrino 
telescopes. 

2 Hawking radiation from PBHs and emis¬ 
sion models 

ffawking showed that a black hole emits particles thermally with a temper¬ 
ature that depends only on its mass, charge and angular momentum. We 
can safely assume that PBffs, formed in the early Universe, are neutral and 
nonrotating because any nonvanishing initial charge or angular momentum 
would have been lost soon after their formation. A black hole emits particles 
with energy in the range (E, E + cLE) at a rate [jl] 

d 2 N _ 1 

dtdE 2nh exp (8nGME/hc 3 ) — (—l) 2s 


T s{E,M) 
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per particle degree of freedom. Here M is the mass of the hole, s is the 
particle spin and T S (E,M) is the absorption coefficient. 

Although Eq. m does not represent a perfect black body spectrum be¬ 
cause T S (E,M) is energy dependent, one can define a temperature 


kT 


he 3 

8v tGM' 


( 2 ) 


The absorption coefficient T S (E,M ) depends on E, M and s and any other 
internal degrees of freedom as well as on the rest mass of the emitted par¬ 
ticle. For relativistic or massless particles F S (E,M) is a function of ME ; 
its exact form has to be calculated numerically [|T0[ . It can be adequately 
approximated in terms of the absorption cross section a s (E,M ) ||: 


r s(E,M) = 


E 2 a s {E,M) 

7T?l 2 C 2 


where o s (E , M) can be replaced by its averaged value 
JdEa s (E,M) 


o* = 


fdE 


56.7 G 2 M 2 /c 4 s = I 
= { 20.4 G 2 M 2 /c 4 s = 1 
2.33 G 2 M 2 /c 4 s = 2 . 


( 3 ) 


( 4 ) 


Particle physics determines the number of states the black hole emits at a 
given temperature. They can be emitted directly by the Hawking mechanism, 
or by the decay and fragmentation products of such particles. The first flux 
we will refer to as “direct”, the latter as “indirect” fluxes. As a consequence 
of emission the black hole loses mass at a rate 

d M a(M) 

= (5) 

where a{M) counts the degrees of freedom of the emitted particles. As 
the black hole radiates, its temperature rises at an increasing rate because 
a(M) increases smoothly at the rest mass threshold for each new massive 
particle ||. For the Standard Model a(M) is shown in Fig. 1 of Ref. 0. The 
details of the final stage of the PBH evaporation thus depend on the high 
energy value of a(M) which is determined by the underlying particle physics. 
Whatever the model, the black hole loses mass approximately as M~ 3 , and 


4 







its temperature increases accordingly. Its luminosity increases roughly as 
M~ 2 as M decreases. We will refer to this runaway process as explosion. 

To a good approximation the solution of Eq. (]3|) for the evolution with 
time of the mass M(t) of a black hole with initial mass M, is 

M(t)^(M/-3 ai ) 1/3 = M.((M) S _y‘ /3 , (6) 

where a = a(M*) and to is the a g e of the Universe. We recall that M* is the 
mass of a black hole, formed in the early Universe and evaporating today. 
One can use Eq. (()|) to relate the initial mass M* of the black hole to the 
lifetime r over which it completely evaporates. 

The most conservative choice for the underlying particle physics model 
is the Standard Model. It provides a lower bound on the high energy value 
of a(M), not only because its degrees of freedom are restricted to known 
particles with the exception of the Higgs, but also because above the decon- 
finement temperature, T D ~ A q C d — 100 — 300 GeV , quarks and gluons 
are emitted as fundamental particles instead of pions and heavier hadrons. 
Quarks and gluons subsequently fragment into (mostly) pions. Only pions 
can be emitted directly below the deconknement temperature, assuming that 
their rest mass is somewhat smaller than A qcd- Other models may include 
more degrees of freedom. For example supersymmetry increases a(M) by 
at least a factor of 3. Instead of trying to list all the possibilities, we only 
mention the other extreme example, the Hagedorn picture, where the num¬ 
ber of degrees of freedom increases exponentially. In Hagedorn-type models 
spectacular explosions are obtained because of these exponentially increas¬ 
ing hadron degrees of freedom. We consider these models unrealistic in the 
framework of QCD and will not consider them further. Throughout this 
paper we will assume the degrees of freedom to be those of the Standard 
Model, even at the highest temperatures, and in our calculations we will use 
T d = 200 GeV, mtop = 170 GeV and 7nm ggs = 100 GeV. 

3 Neutrino fluxes 

We are now ready to calculate the neutrino fluxes for each species. The direct 
neutrino flux, which is the same for all neutrino flavors, is given by Eq. ©• 
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The indirect neutrino fluxes are flavor dependent and fall into two categories. 
They can be the decay products of directly emitted heavier particles, like 
muons or taus. Above the QCD deconfinement temperature they can also be 
produced by the decay of quark and gluon fragmentation products, mostly 
pions in quark jets. These fluxes are studied in detail in |J. Here we shall 
make the adequate approximation that all the fragmentation products are 
pions, hence the decay of fragmentation products only contributes to the 
muon and electron neutrino flux through charged pion decay. 

Our main results are independent of the explicit choice of the quark frag¬ 
mentation function. We use the empirical fragmentation function 0 


dAL 15 


d^ 


= —(z-i ) 2 z~ 3/2 


16 


(7) 


where z = E v /E and E is the quark (jet) energy. This choice of the fragmen¬ 
tation function is convenient because it yields realistic, analytic expressions 
for the hadron decay fluxes when we further approximate the flux by a delta 
function at its peak value 0. 

Two thirds of the pions are charged and decay into u tl +// and their charge 
conjugates. The direct quark emission is convoluted with the fragmentation 
function and the pion decay distributions. Naturally, whenever fragmenta¬ 
tion is considered, the flux depends on the PBH mass at the moment of 
emission. The muons from pion decay also decay (/i _ —> v e e~), and thus 
contribute to the and the v e fluxes. To calculate the contributions of such 
decays the convolution must therefore also involve muon decay. 

The decay distribution of a particle from any decay is given by: 


dn 1 dT 

d E = Fd E ’ 


where T is the decay width. The distribution is normalized to the number of 
identical particles produced by a single decay. The calculations are straight¬ 
forward. Below we explicitly list the distributions used because they correct 
errors in an earlier calculation [§ : 

Pion decay: —> / u ± ^(w)- The z/ /x distribution is given by: 



d E v 



™l ~ m l ) \ 

2m, V ’ 


(9) 
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1 


( 10 ) 


in the rest frame of the pion (/3 = 0 ), and 


dn^ _ ml 
d E v ml - ml j £2 - ml ’ 


where the allowed range of E v is 

| {ml - ml) < F < \ {ml - m 2 ) 
E, + - ml ~ E„- 


( 11 ) 


in the lab frame (/? ^ 0 ). 

Muon decay: ^ u e (u e ) e ± The distribution of the muon neu¬ 

trinos is given by: 


d n v 


I 5 _ 3 2 , 2 3 

/? L 6 2 C ^3 


d E v 7 


x 


2 (1 + /3) 2 [3e 2 -|(^f) e ; 

Similarly the electron neutrino distribution is: 

d [1 - 3e 2 + 2e 3 ] 


d n v 


x 


^ [ 4(l+/3) 2 [3e 2 -(^f ) e 3 

In both the above cases e = E u /E™ ax , where 


= —pf(l +/?), 


for l=f < £ < 1 , 


for o < £ < y 


for j+§ < £ < 1 , 


for 0 < £ < y 


( 12 ) 


(13) 


(14) 


is the maximum neutrino energy, and the electron mass is neglected. 

Tau neutrinos represent a somewhat different story. They can hardly 
be produced in hadron decays since decays of hadrons into v T are highly 
suppressed. Decays into v T V T proceed via Z° exchange and therefore are 
suppressed with respect to decays mediated by a photon. Decays into a 
r — u T or r — f pairs are kinematically forbidden except for heavy hadrons 
which rarely appear as quark or gluon fragmentation products. Moreover, 
the decays into r — v T channels have small branching ratios. In the end most 
of the v T flux is direct. The only other significant contribution comes from 
the decay of direct r’s. 


7 



For the neutrino fluxes from r decay we estimate that ~ 40% of the 
decays occur in leptonic modes, with the rest in semi-hadronic modes. In 
the free quark approximation this comes from the relative number of colors 
(3) to neutrino flavors (2) into which the decay may proceed. The number 
distributions for each type of particle are the same as the above three body 
distribution in the limit of massless final states. Thus an approximate v T flux 
is given by replacing m M by m T in the above formulas. The v e and z/ /; fluxes 
coming from r —> ev e v T and r —> respectively, are each 20 % of the 

total v T flux from r decay. 

In the following sections we will compute f(x, M), the total neutrino flux 
of energy E = ( 8nGM/hc 3 ) x for a given species. For z/ M , the direct, the 
muon decay, and the fragmentation flux contribute to f(x,M). The direct 
flux depends on x only, while the fragmentation flux depends also on the 
PBH mass M at the time of the emission. 


4 Diffuse neutrino fluxes 

Our calculation of the diffuse neutrino spectrum closely parallels that of 
Ref. pi for gamma rays. In order to calculate the present neutrino flux of 
energy Eq we must integrate the contributions over cosmological time t for 
all PBH of mass M at the blue shifted energy E 

E=—, (15) 

r 

where r = R/Ro = (1 + Z ) -1 , Z is the redshift at time t. The subscript 0 
denotes the value of a quantity at the present epoch. We assume a standard 
Friedman cosmological model. Thus the Universe is radiation dominated for 
t < 10~ 6 s and matter dominated for t > 1CU 6 s. 

Assuming a uniform distribution of PBHs, created shortly after t — 0, 
the initial number of PBHs per comoving volume of mass less than M, is 

/■Mi/M, 

n(Mi) =M / d y s(y) , (16) 

Jo 

where y = M//M* and s(y) is a dimensionless function that determines the 
form of the initial mass distribution. Here s(l) = 1. J\f is the initial number 
density of PBH per logarithmic mass interval at M — M*. From various 


cosmological considerations we know that Af < 10 4 pc -3 . PBH clustering 
in the galaxy can increase their local density by as much as a factor of 
£ = 10' which dramatically improves our chances to observe them [9]. In our 
calculations we will not assume any clustering, ie. £ = 1. 

If PBHs form from scale invariant density perturbations with a Gaussian 
spectrum, then 

s(y) = y~ p , (17) 


where 2 < (3 < 3, depending on the equation of state of the Universe at the 
formation of PBHs []3[. If p — yp, where p is the presure and p is the energy 
density, then 


P 


19-37 
1 + 7 


+ 1 . 


(18) 


Let f(x, M ) represent the total neutrino flux of energy E u as defined in 
the previous section, then the diffuse flux per unit area today is: 


d 4 N Afc r j o r 

dA d£l dt dE M = L ^ + Z) j Ay s(y) , (19) 

where y is to be integrated over all values of the initial masses of the PBHs 
present at t', and x is the value of (87 t GM/hc 3 ) E u at t! and y. Absorption on 
ambient matter in the Universe should be taken into account. For neutrinos 
this is negligible. Making use of (|J) and (!17j), the integral ([T9]) becomes: 


d A N v 

dA dQ dt dE u0 


Afc 
47nr 3 


dt’ 



dx x 2 


xf(x,M) 




-(£± 2 ) 

5 


( 20 ) 


where x* = (87T GM*/hc 3 ) E u . The right hand side can be integrated numer¬ 
ically once the contributions to f(x, M) are specified. For electron neutri¬ 
nos we include the contribution from the decay of secondary muons (from 
the decay of fragmentation pions), the decay of direct muons and taus and 
the direct flux. For muon neutrinos we include the contributions from the 
decay of fragmentation pions, the decay of secondary muons, the decay of 
direct muons and taus and the direct flux. In the fragmentation case the 
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x-integration limits are modified. For tau neutrinos only the direct flux and 
the flux from direct tau decay contribute significantly as discussed earlier. 

In Figs. 1 and 2 we plot the resulting fluxes. The fluxes have been cal¬ 
culated for a PBH density taken to be the maximum allowed by bounds 
derived from gamma ray measurements 0. The main characteristics of the 
neutrino fluxes are: i) in the low energy region the slopes of the curves are 
/3 dependent, while at high energies they decrease as FV 3 , independent of (3 , 
ii) all fluxes display a “knee” (a change of slope) in the 50-100 MeV region 
corresponding to the peak emission energy of black holes of mass M*, iii) 
the fragmentation spectra display a knee at around 10 MeV and, finally, iv) 
the direct spectra are negligible above the knee and dominate in the region 
below. 


5 Background fluxes and bounds on PBH 
abundances 

It is in principle possible to obtain an improved bound whenever the pre¬ 
dicted neutrino fluxes, computed for the maximal PBH abundance allowed 
by the gamma-ray bound, exceed the neutrino backgrounds. These are pre¬ 
dominantly of solar and atmospheric origin and can be separated by energy 
region. 

The solar flux extends to 18.8 MeV where it is sharply cut off. Below the 
cutoff it dominates with a very featured spectrum which traces the different 
nuclear reactions in the Sun’s interior. The reaction yielding the lowest flux 
but the most energetic neutrinos is He + p —> v e + ..., usually abbreviated 
hep. This flux is large compared to the diffuse PBH flux but it has two 
distinguishing characteristics that can be used to suppress it: i) it is of u e 
flavor and ii) directional. The required rejection rate of order 10~ 6 — 10 -9 
is challenging in spite of the excellent sensitivity of the detectors to the 
neutrino direction, especially for the hep and ' B fluxes. Flavor identification 
at this level seems difficult since ry or v T are not detected. Finally, flavor 
identification is only valid in the absence of solar flavor oscillations, an issue 
which is unsettled. 

Above the solar cutoff the background neutrinos fluxes are dominated by 
the decays of hadrons produced by interactions of cosmic rays with atmo- 
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spheric nuclei. The main contribution is from pions with v e : ~ 1 : 2. The 

atmospheric neutrino fluxes are understood with roughly 10% precision. At 
the lower energies they are strongly dependent on the magnetic latitude and 
arrival directions because of the geomagnetic cutoff for cosmic rays. Loca¬ 
tions with low geomagnetic latitude have smaller atmospheric neutrino fluxes 
at low energies which is an advantage if a bound on PBHs is to be extracted. 

The atmospheric flux of v T is multiply suppressed. The open decay chan¬ 
nels involve heavy parents such as D s or bottom hadrons JT2| which have 


smaller production cross section than pions and kaons and branching ratios 
for v t t decays which are over an order of magnitude (4% for D s ) below those 
for pions and kaons. Production of heavy hadrons has a higher threshold 
energy. Because of the steepness of the cosmic ray parent flux, increased 
thresholds result in a drastic suppression of the flux at low energies. The 
calculation is very similar to that of the prompt muon neutrino flux of charm 
origin [|13j and subject to similarly large uncertainties associated with our 
poor understanding of the production cross section for D s in the region near 
threshold. There is no point then in performing a detailed calculation, and 
we obtain a simple estimate of the v T flux by rescaling the prompt muon 
flux by a factor. Considering that D s production is suppressed and that the 
branching ratios of the decays D s —> u T +... are smaller than the correspond¬ 
ing ones for the decays of other Ds, we approximate the above factor to be 


of order 0.1. This agrees with the estimates made for LHC [|14[ . 

In Fig. 3 we compare PBH fluxes to the atmospheric u e and fluxes at 
a north U.S. location [21], to the solar neutrino hep flux and to our estimate 
of the atmospheric u T flux. The PBH fluxes are calculated for [3 = 3, because 
one can assume that at the formation of PBHs the Universe is radiation 
dominated. The figure allows one to identify the most promising regions for 
obtaining an improved bound. The high energy region is promising since 
the background falls faster than the signal. Unfortunately prompt neutrinos 
from charm decay are expected to dominate at the highest energies over the 
fluxes of pion and kaon decay origin and may thus prevent the atmospheric 
neutrino flux to fall below the PBH flux. On a more speculative note, the 
high energy neutrino flux may be dominated by point sources which is not a 
problem unless a large number conspire to produce a diffuse flux. Searching 
for a v T flux seems to be most promising. Even nonobservation is of interest 
as it may yield important bounds on PBH abundances. 


11 












6 Final stage evolution: detection of z/s from 
individual black holes 


As we have seen in Sect. 2 the final stage of a PBH is explosive. In order to 
determine the present density of these exploding black holes one must look 
for individual sources. In this section we show that high energy neutrino 
telescopes can significantly contribute to this search. Their threshold energy 
is of order E > 10 GeV for neutrinos. The main contribution to the flux 
f(x,M ) comes, in the case of explosions, from the decay of fragmentation 
pions. This is because other contributions are either smaller or peak at lower 
energies. 

Because we are going to consider PBHs whose lifetimes are short com¬ 
pared to a typical observation time (of order of years), the relevant quantity 
is the time integrated flux 


d N v 
d E 



d tf{x,M) 


1.19 10 ^GeV- 1 


n GeV 
V E 


3 



d xx 2 a 1 f(x,M), 


( 21 ) 


where x = E/kT. Here we have used Eqs. (|2|) and ([5]) to convert the time 
integral into an integral over the inverse temperature. We can neglect redshift 
effects because we will find that one does not expect to see PBH explosions 
from further than a few parsecs. 

The interact with the Earth and the secondary upcoming muons pro¬ 
duce the experimental signature in the underground detector. Only muons 
entering the detector from below are detected, because the background of 
downward atmospheric muons dominates any neutrino-induced signal. The 
probability that a neutrino produces a muon above the threshold energy E D 


at the detector 115 


is 


Pu^ti(E u , E d ) — N a f 

J E 


E " dE^iEp, E U )R,(E„ E d ), (22) 

Ed 


where N A is Avogadro’s number, da/dE^ is the differential charged current 

and R M is the muon range fH|. We updated da/dE^ 


cross section [|1S 

of 0, [L9| by using the latest small-a; parton distribution functions of 
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The evaluation of Eq. 
routine VEGAS ||T7|. 


was carried out using the Monte Carlo integration 


The predictions relevant to experiment are obtained by integrating the 
muon flux in the detector above E n . The number of observed events is given 
by: 


A c°° dN 

N u (E d , 9) = — 2 d E P^exp [~a tot (E u )N A X(9)} —A (23) 

dwr z Je d d E 

where r is the distance to the PBH, A is the sensitive area of the detector, the 
exponential factor accounts for absorption of neutrinos along the chord of the 
Earth X{9) and <7tot is the total cross section |19|. Here 6 is the zenith angle 
(0 < 6 < 90°). The dominant contribution to the integral (^3j) comes from 
and is evaluated for nearly vertical neutrinos (6 ~ 0), because only neutrinos 
with small 9 travel far enough in the Earth to interact and produce a muon. 
Detectors under construction may, eventually, reach an effective area of order 
A ~ 1 km 2 . The threshold energy can be anywhere between a few GeV to a 
few TeV. 

Because the number of events is inversely proportional to r 2 , there is a 
maximum r for which the PBH is detectable. This distance depends on the 
initial temperature of the PBH and is increased for lower temperatures be¬ 
cause of the increased observation time. Fig. 4 shows this maximum distance 
as a function of the detector threshold energy E D for initial temperatures 
ranging from 50 GeV to 20 TeV. The limiting distance is determined by the 
requirement that we detect at least 10 muons per burst. 

For a PBH to be detectable it has to stand out over background. In 
this case the high energy atmospheric flux is the relevant background. Our 
calculation follows Ref. p2|, and is in agreement with other calculations and 
measurements of the atmospheric neutrino background pi 
background events in the detector is: 


The number of 


IVbgr = r A 


I Ed 


dE ME) dA dt ^ dE , 


where iv D is the angular resolution of the detector 

0.314 sr 


uj d (E) = 


E{ GeV) ’ 


(24) 


(25) 
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and r is the lifetime of the black hole. For detectability we require that the 
signal ([ 


) exceeds the background (24) by 5cr. The maximum distance for 
which this condition holds, is shown on Fig. 4 as a function of the detector 
energy for initial temperatures ranging from 50 GeV to 20 TeV. 

The distance of an observable PBH must fall below both curves on Fig. 4. 
For kTi > 1 TeV the 10 muon signal condition is stronger. This is because the 
lifetime of such a PBH is so short that few background events can accumu¬ 
late during its observation. For lower initial temperatures the 10 muon signal 
condition is stronger for high detector thresholds, while the 5cr condition is 
stronger for low thresholds. This combined condition is shown on Fig. 5. The 
optimal condition is obtained for a PBH with kTi — 100 GeV. The lifetime 
of such a PBH is about 1.5 months, longer than the lifetime of PBHs usually 
referred as “explosive” in the literature, but still small compared to the de¬ 
tection time. The maximum distance is 2.3 x 10 _3 pc at E D ~ 140 GeV, which 
gives an upper bound 1.2 x 10'pc _ 3 yr ^ 1 on the number of PBHs exploding 
in our neighborhood per year. Here we used detector solid acceptance angle 
27r, observational time 3.3 years, and assumed a uniform PBH distribution 
in our neighborhood. The detector acceptance solid angle is 2n because only 
upward muons have experimental signature. In order to compare our result 
to existing experimental bounds, we choose the observation time to be that 
of the search establishing the strongest bound. It is 8.5 x 10 5 pc - 3 yr _1 from 
high energy 7 ray measurement from kTi — 10 TeV PBH, and it is of two 
order of magnitude smaller than that of any other direct search 
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Our result is stronger than any existing experimental bounds, except 
that of |23|]. Our bound is about an order of magnitude larger, but can 


be interesting if for some reason detection of neutrinos is preferable over 
detection of 7 rays from PBHs. The issue has indeed been raised whether 7 
rays are emitted as independent particles in the presence of thermal electrons 
surrounding the black hole |24] . 


7 Detector considerations 

It obviously remains to be seen whether detectors will ever reach the sensitivi¬ 
ties to observe the neutrino fluxes predicted in our calculations. In the diffuse 
case, the challenge is formidable as can be seen from Fig. 3. One must im¬ 
prove on the sensitivity to detect atmospheric neutrinos by over three orders 
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of magnitude even when restricting observations to the favorable vertical neu¬ 
trinos. The most promising search involves tau neutrinos with energy below 
10 GeV. Detecting the r-lepton emerging from a charged current interaction 
may be possible since the neutrino energy can exceed the r-lepton mass but 
this will certainly represent a challenge. Alternatively one may look for the 
difference between the total neutrino flux and those of electron and muon 
flavor. This procedure requires a precision of 0.1% which is probably out of 
question in the near future. The search for explosions presents us, in the 
end, with the most realistic opportunities. 

8 Conclusion 

In this paper we have calculated the neutrino fluxes from PBHs and their 
detectability for both low and high energies. We have updated the previous 
calculations by including contributions not only from direct neutrino emis¬ 
sion, but from muon and tau decays and from quark fragmentation. The 
latter is especially important because of the large number of quark degrees 
of freedom. 

In the diffuse case the fluxes are increasing functions of (3 in the allowed 
range 2 < /3 < 3. The /y, and v e fluxes are dominated by the fragmentation 
contribution for E v less than about 100 MeV. Around 1 MeV it is about an 
order of magnitude higher than the other contributions. For larger energies 
all contributions decrease as E~ 3 , and the total flux is dominated by the 
direct emission of neutrinos. In the case of u T the direct flux dominates at 
all energies, because there is no relevant fragmentation contribution. In the 
individual PBH case the time integrated flux is relevant. The contribution 
from the first decay of fragmentation pions is the largest. 

We have compared the calculated fluxes to the relevant background to 
determine which neutrino species at what energy range has large enough flux 
to be detectable. In the diffuse case the v e and v M fluxes are larger than the 
atmospheric background only for E v smaller than a few MeV. For v e the huge 
directional solar hep background can cause further observational difficulties. 
Most promising is the z/ r case, because the only relevant background, the 
atmospheric z/ T flux, is highly suppressed. 

Our best chance is to look for individual PBHs in our neighborhood. 
We have shown that the new kilometer-scale, high energy ( E D ~ 1 GeV — 


15 



1 TeV) neutrino telescopes can play a significant role in this search. They are 
sensitive to v^s. Furthermore at these high energies the only background, 
the atmospheric neutrino background, is suppressed. The bound that can 
be established with these detectors is comparable to existing experimental 
bounds, obtained by searching for high energy 7 ray bursts. 
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Figure captions 


Fig. 1. Contributions to diffuse neutrino flux (a) for z/ e , (b) for z/ M and (c) 
for v T . 

Fig. 2. Total diffuse neutrino fluxes for (3— 2, 2.5 and 3 and (a) for z/ e , (b) 
for and (c) for u T . 


Fig. 3. Total diffuse neutrino fluxes for (3 —3 compared to the atmospheric 
neutrino backgrounds and to the solar hep flux. 


Fig. 4. The 10 muon signal (dotted lines) and 5 a signal (solid lines) con¬ 
ditions as the function of the detector threshold energy, E D . Dotted lines 
correspond to kTi = 50 GeV, 100 GeV, 200 GeV, 500 GeV, 1 TeV, 2 TeV, 
5 TeV, 10 TeV and 20 TeV from top to bottom, solid lines correspond to the 
same kTi values from bottom to top. 


Fig. 5. Combined 10 muon signal and 5 a signal conditions as the function 
of the detector threshold energy E D for kTj = 50 GeV, 100 GeV, 150 GeV, 
200 GeV, 500 GeV, 1 TeV, 2 TeV, 5 TeV, 10 TeV and 20 TeV. 
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